Variational Monte Carlo Study of Anderson Localization in the Hubbard Model 
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We have studied the effects of interactions on persistent currents in half-filled and quarter-filled 
Hubbard models with weak and intermediate strength disorder. Calculations are performed using a 
variational Gutzwiller ansatz that describes short range correlations near the Mott transition. We 
apply an Aharonov-Bohm magnetic flux, which generates a persistent current that can be related 
to the Thouless conductance. The magnitude of the current depends on both the strength of the 
screened disorder potential and the strength of electron-electron correlations, and the Anderson 
localization length can be extracted from the scaling of the current with system size. At half 
filling, the persistent current is reduced by strong correlations when the interaction strength is large. 
Surprisingly, we find that the disorder potential is strongly screened in the large interaction limit, 
so that the localization length grows with increasing interaction strength even as the magnitude of 
the current is suppressed. This supports earlier dynamical mean field theory predictions that the 
elastic scattering rate is suppressed near the Mott transition. 



I. INTRODUCTION 

The Hubbard model is the standard model for 
strongly-correlated particles in a tight-binding lattice, 
and has been widely useful as a minimal model to de- 
scribe a variety of transition metal oxidesji quantum 
dot arrays^ and, more recently, trapped atomic gases. 3 
Real materials are disordered and extensions of the Hub- 
bard model, most notably the Anderson-Hubbard model 
(AHM), have been introduced to study the subtle inter- 
play between disorder and interactions. Like the Hub- 
bard model, the AHM has not been solved exactly, ex- 
cept in special cases, and appears to have a rich phase 
diagram. 

Over the past decade, there has been considerable 
interest in Anderson localization in strongly-correlated 
metals near half filling. At half filling, Coulomb interac- 
tions can induce a gapped "Mott insulating" state, and 
the question is how states near the Mott state are af- 
fected by strong interactions. This is particularly rele- 
vant to the cuprate high temperature superconductors, 
where the electronic properties are tuned by chemical 
doping. Chemical doping, in addition to changing the 
carrier density (and therefore proximity to the Mott 
state), disorders the materials. Experiments have gener- 
ically found states exhibiting characteristics of Anderson 
localization that lie between the metallic and Mott in- 
sulating phases The location of the metal insulator 
transition (MIT) appears to be a function of both the 
itinerant hole density and the level of disorder, but not 
the residual (T — > 0) resistivity, 4 which hints at a possible 
failure of the conventional one-parameter scaling theory^ 
More generally, these experiments raise the possibility 
that the Anderson MIT in strongly correlated materials 
may differ qualitatively from that in conventional weakly 
correlated metals. 

In this work, we use the AHM to study the effects 
of strong correlations on Anderson localization. The 



Anderson-Hubbard Hamiltonian is 



H 
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+ (UhifUii + Vjhj) , (1) 



where the indices i and j refer to sites on a d-dimensional 
lattice, where c\ a annihilates a spin-cr electron from site i, 
hi a = c^ ia Ci a is the number operator for site i, and where 
hi = yVniq. We restrict ourselves to nearest- neighbor 
hopping, so that = — t if i and j are nearest neighbors 



and tj 



otherwise. Disorder is introduced via the site 



energies Vi, which are chosen randomly from a uniform 
distribution of width W (so -W/2 < V t < W/2). The 
AHM is therefore characterized by three energy scales: 
the kinetic energy t, the short range Coulomb interaction 
U, and the disorder strength W. 

A rough schematic phase diagram for the half-filled 
AHM in three dimensions is shown in Fig. [T] (the phase 
diagram in lower dimensions has not, to our knowledge, 
been established). This phase diagram is based on a 
number of published calculations^— which, while they 
disagree on details, agree on general features. Note that 
we have chosen to show the simplest nonmagnetic phase 
diagram (in which magnetic phases are suppressed) and 
have omitted phases that are not universally seen. 

There appears to be consensus on the strong disorder 
limit of the phase diagram: as a function of increasing 
U, there is a continuous transition from a gapless An- 
derson insulator to a gapped Mott insulating state at 
U = U C (W) ~ W An interesting, recently-discovered 
wrinkle in this picture is the existence of a kinetic energy- 
driven zero bias anomaly in the density of states^ - — 
which was found in the one and two dimensional AHM 
and should presumably be present in higher dimensions. 

There is less consensus on the phase diagram at weak 
and intermediate disorder. Some calculations have pre- 
dicted that the density of states vanishes continuously at 
the MIT in two 1 ^ and thre e) 10 ' 19 dimensions while dynam- 
ical mean field theory 9 - 2 ^ (DMFT) finds that the den- 
sity of states is discontinuous across the MIT, and varia- 
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FIG. 1: Schematic phase diagram for the Anderson-Hubbard 
model in three dimensions at half filling. Thin lines indicate 
continuous transitions between metallic, Anderson insulating 
and Mott insulating phases; the thick solid line indicates a 
region of the phase boundary where DMFT predicts a dis- 
continuity in the density of states. 



tional Monte Carl o 22 ' 23 calculations predict a first order 
Mott transition. Although experiments are often com- 
plicated by the occurrence of broken-symmetry phases, 
paramagnetic Mott transitions in real materials are gen- 
erally found to be consistent with DMFT>2i The details 
of the T = MIT are important since they affect the 
electronic properties of the metal adjacent to the Mott 
transition. 

To illustrate this point, we consider a simple weak- 
coupling expression for the screened random potential 
Vi. Within mean-field theory, Vi = Vi + \U{m — no), 
where hq is the spatially-averaged charge density and 
rii = (hi). For weak disorder, rij = n +J2j Vjdrii/dVj. A 

rough estimate for the charge susceptibility is drii / dVj ~ 
—2pQ§ij, where po is the single-spin density of states at 
the Fermi energy. Then 



Vi 



V 



1 + PqU 



(2) 



In a weak-coupling metal, po depends only weakly on 
U and the overall effect of increasing U is to screen the 
disorder potential. On the other hand, in a gapped insu- 
lator, po — an d the impurity potential is unscreened. 
While Eq. (J2)) applies to weakly-correlated insulators, we 
also expect screening to vanish in (strongly-correlated) 
Mott insulators because the charge is unable to rearrange 
itself in response to the disorder potential. The screened 
potential should therefore be a nonmonotonic function of 
U , obtaining a maximum value somewhere between small 
U (where Eq. ([2} holds) and the critical value U C (W). 
Indeed, several calculations have found that the impu- 
rity potential becomes unscreened as U — V U c for strong 
disorder 

On the other hand, DMFT predicts that, for weak dis- 
order, a narrow quasiparticle band develops at the Fermi 



energy near U c , such that the density of states remains 
finite up to U c . It has been shown that because this 
band tends to be pinned to the Fermi energy, impurity 
scattering within the band is reduced i 26 ' 27 In fact, for a 
particle-hole symmetric band, a weak impurity potential 
is predicted to be perfectly screened as U —> U c . This re- 
markable result has, to our knowledge, only been found 
in DMFT calculations and has not been confirmed by 
other theoretical techniques. 

In this work, we study persistent currents in an 
Anderson-Hubbard lattice threaded by an Aharonov- 
Bohm flux and focus on the region of strong correlations 
near the Mott transition. These currents are affected by 
both the disorder potential (which tends to localize the 
quasiparticles), and by strong correlations (which reduce 
the quasiparticle spectral weight). One of the main ob- 
jectives of our analysis is to separate these two competing 
effects, and we find the paradoxical result that the local- 
ization length may grow with increasing U, even as the 
magnitude of the current is suppressed. 

In Sec. HU we describe how we calculate the persistent 
current for a ring or torus threaded by an Aharonov- 
Bohm flux. In Sec. IIIII we compare clean and disordered 
systems with and without strong correlations, from which 
we develop a qualitative sense of the relative importance 
of disorder screening and strong correlations. We then 
show how these affect the localization length, which is 
extracted from the scaling of the persistent current with 
system size. We present results for two cases: quarter 
filling where correlations are relatively weak, and half 
filling where correlations are strongest. We find that, 
while the quarter-filling results are consistent with one- 
parameter scaling, there is evidence for a breakdown of 
one-parameter scaling at half-filling when the interac- 
tion strength U is greater than the disorder strength W . 
These results are discussed in Sec. IIV1 



II. CALCULATIONS 

The Hamiltonian is given by Eq. ([1]). An Aharonov- 
Bohm flux <I> is introduced through a complex phase in 
the hopping matrix elements 



U 



-t exp 



il-KOLXi 



(3) 



where < a < 1, L is the circumference of the ring (in 
one dimension) or torus (in two dimensions), and Xij — 
±1 is the electron displacement in the x direction during 
the hop from site j to site i. The parameter a is the 
magnitude of the Aharonov-Bohm flux in units of the flux 
quantum «3>o = hc/e. The magnitude t of the hopping 
matrix element between nearest-neighbor sites is taken 
to be the unit of energy throughout this work (ie. t = 1). 

The variational ground state wavefunction has the 
Gutzwiller form 



I* 



gwf) = fblVVs) 



(4) 
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FIG. 2: Variational energies of the PMGW and DFSGW 
states for a one dimensional chain with L = 50 and W — 4t. 
Results are at half filling, and are for a single configuration 
of disorder. The energy per electron is shown. Inset: En- 
ergy difference AE = Sdfsgw — £pmgw between the two 
variational states. 



where |i/> ps ) is a Slater product state wavefunction and 
Pg = njl — (1 — flO^if^il.] is the Gutzwiller projection 
operator. The parameter g is a variational parameter 
that introduces correlations into the wavefunction; the 
limit g — corresponds to a state with no double oc- 
cupancies, while g — 1 gives an uncorrelated wavefunc- 
tion I^gwf) = hAps)- This simple Gutzwiller wavefunc- 
tion cannot, without additional Jastrow factors, describe 
the Mott transition) 28 ' 29 nonetheless, the Gutzwiller pro- 
jection is useful as a qualitative tool for understanding 
strong correlation physics near the Mott transition. 

It is usual to take |"0ps) to be the many-body ground 
state of some effectively noninteracting Hamiltonian, and 
we consider here two cases: (i) l^ps) is the self-consistent 
paramagnetic Hartree-Fock ground state of H, and (ii) 
\ip ps ) is the ground state of a noninteracting Hamiltonian 
with a screened disorder potential, 



H* = E 1 



(5) 



The parameter e is treated as a variational parame- 
ter, on the same footing as g. We call the projected 
wavefunctions obtained from (i) and (ii) the paramag- 
netic Gutzwiller wavefunction (PMGW) and the disor- 
dered filled sea Gutzwiller wavefunction (DFSGW) re- 
spectively. 

While it is possible to include more variational pa- 
rameters, for example by taking g and e to be site 
dependent,—^ we found previously^ that it is sufficient 
to include spatial inhomogeneity only in \ip ps ) when the 
disorder is not too strong; it is only for strong disorder 
(relative to the bandwidth) that a spatially varying g and 
e are essential. This represents a large computational 
savings. 



The parameters g and e are determined by minimizing 
the energy functional 



GWF 



HI* 



GWF 



GWF I* GWF 



(6) 



For large systems, £ cannot be evaluated exactly, and 
a variational Monte Carlo (VMC) method is used. The 
idea is to expand \i[>ps) in terms of Fock states \n) (states 
with a definite number of electrons on each site), such 
that \ipps) = En a «l n )- It is then simple to apply the 
Gutzwiller projection: PgIV-'ps) = En^™! 71 ) wnere <5 m = 
g Dm a m , with D m the number of doubly occupied sites in 
the ket 



\m). 



£ 



Then the variational energy can be written 
= Em I 5 "*! 2 En (H mn a n /a m ) 

E P M 2 



The double sum in the numerator is computationally pro- 
hibitive to calculate and, in the simplest approach, the 
sums over m and p are performed approximately using 
the Metropolis algorithm with weighting factors \a m \ 2 
and \a p \ 2 respectively. More sophisticated VMC algo- 
rithms exist, and in this work, we used the modified 
Metropolis algorithm described in Ref. HH. 

The results of the energy minimization are shown in 
Fig. [2l and suggest that there is little difference between 
the PMGW and DFSGW states. We shall see in the next 
section that although these two product states produce 
almost identical results for small L, they make quanti- 
tatively different (although qualitatively similar) predic- 
tions for the persistent current when L is large. This 
difference stems from how screening is handled: in the 
PMGW approximation, the self-consistent charge densi- 
ties are obtained from a Hartree-Fock calculation before 
£ is minimized with respect to g; in the DFSGW approx- 
imation, £ is minimized with respect to the projection g 
and screening e simultaneously. We expect, therefore, 
that the PMGW state provides an upper bound on the 
screening, and where we have been able to compare with 
published results the DFSGW gives a more accurate es- 
timate of the persistent current. 

The Aharonov-Bohm flux in the Hamiltonian generates 
a persistent current 



J = 



GWF 



GWF 



(^gwfI^gwf) 



with 



itea 



E 

(i,j),o 



Xij exp 



C ia C jcr- 



(7) 



(8) 



We note that the persistent current is approximately re- 
lated to the Thouless conductance via gxh ~ J/ot- 

The persistent current depends on the length L of the 
chain or torus, on the interaction strength, and on disor- 
der. In the next section, much of the discussion involves 
a comparison of J for systems with and without disorder. 
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FIG. 3: Distribution of the persistent currents for an ensemble 
of one dimensional chains with L = 50, W = it and U = 8t. 
Here and throughout this work, the Aharanov-Bohm flux is 
a — 0.25. Results are shown for paramagnetic Gutzwiller 
(left) and paramagnetic Hartree Fock (right) approximations. 
Smooth curves are Gaussian fits to the data for comparison. 



However, for one dimensional rings, it is also possible to 
perform a finite size scaling analysis from which the local- 
ization length can be extracted. For localized electrons, 
the persistent current satisfies J = Joexp(— L/£) in the 
limit L > where £ is the localization length, and it 
is thus possible to extract the localization length from 
logarithmic plots of the persistent current as a function 
of system size. 

The current J depends on the details of the disorder 
configuration and, in the noninteracting case, the distri- 
bution is log-normal when the system size is much larger 
than the localization length. 32 Examples of the distribu- 
tion of the current values are given in Fig. [3] Because 
of the approximate log-normal current distribution, we 
calculate the localization length from the typical value of 
the current 



Ji 



typ 



exp (in J) 



(9) 



where the bar denotes an average over disorder configu- 
rations. 



III. RESULTS 



A. The Persistent Current 




FIG. 4: Persistent current as a function of U for PMGW and 
DFSGW ground states for both a clean (W = 0) and disor- 
dered (W = 8t) two dimensional plane. Results are shown 
at half filling (n = I) for a single disorder configuration with 
linear dimension L — 8. Results with g set to 1 by hand in- 
clude screening but not strong correlations; these are labelled 
PMHF (paramagnetic Hartree-Fock, obtained from PMGW) 
and DFS (disordered Fermi sea, obtained from DFSGW). The 
DFS curves have the same value of e as the DFSGW curves. 
The inset shows the dependence of the variational parameters 
on U. 




FIG. 5: As in Fig. [4] but for a one dimensional chain with 
L = 50 and W = it. The inset shows PMGW and DFSGW 
results for a short [L = 8) chain. 



We begin by showing results for a single configuration 
of disorder as a function of U. Figures 0] and [5] show the 
persistent current at half filling as a function of the inter- 
action strength in two and one dimensions respectively. 
The disorder strength is chosen to be equal to the band- 
width of the disorder-free lattice, namely W — 8t in two 
dimensions and W — At in one dimension. Results are 
shown for both the PMGW and DFSGW ground states. 
In the two dimensional case, the two approximations pro- 
duce nearly identical currents, while in one dimension 
they are quantitatively different, although qualitatively 



similar. The distinction appears to be that L <C £ in 
Fig. 2] while L 3> £ in Fig. [SJ in support of this, we have 
found that the DFSGW and PMGW currents in one di- 
mension are nearly the same when L is small (inset of 
Fig. p. 

The main point of figures S] and [5] is that they allow us 
to distinguish the effects of screening and strong correla- 
tions on the persistent current. To separate these effects 
we show, in addition to Jpmgw and JdfsgW: results for 
calculations without disorder (W — 0), without strong 



5 



correlations, and with neither disorder nor strong corre- 
lations. 

To suppress strong correlations, the Gutzwiller projec- 
tion is explicitly turned off. These calculations are la- 
belled "paramagnetic Hartree-Fock" (PMHF) and "dis- 
ordered filled sea" (DFS), and are obtained by setting 
.9 = 1 in the PMGW and DFSGW calculations respec- 
tively. (In the case of the DFS calculations, the screening 
e is the same as in the corresponding DFSGW calcula- 
tions.) Figures 0] and [5] show that, as U increases, both 
Jpmhf and Jdfs increase monotonically, and approach 
J of the disorder- free system (PMHF with W = 0). This 
indicates that, in the absence of strong correlations, dis- 
order is progressively screened with increasing U. In 
Fig. El the PMHF curves approach the disorder-free limit 
faster than the DFS curves, indicating that screening is 
stronger in the PMHF calculations; however, the quali- 
tative trends are the same in both cases. 

When the Gutzwiller projections are included, the cur- 
rent is nearly identical to the unprojected current for 
small U, and begins to differ from the unprojected cur- 
rent when U rs W. We have checked several different 
values of W and found that this is generally the case: 
for intermediate disorder strengths, W is the crossover 
scale beyond which strong correlations begin to play a 
significant role. When U > W, the persistent current 
is progressively suppressed as U increases. Within the 
Gutzwiller approximation, this suppression is interpreted 
in terms of a quasiparticle spectral weight that is reduced 
as the Mott transition is approached. This mechanism is 
well-known and has been proposed, for example, to ex- 
plain the anomalously low superfiuid screening currents 
in underdoped cuprate superconductors. 33 

The most remarkable feature of Figs. H] and [5] is that, 
although the Gutzwiller projected currents Jpmgw an d 
Jdfsgw are suppressed by strong correlations at large U, 
they approach the value of the current in the (strongly 
correlated) disorder-free system. This indicates that, in 
spite of the strong correlations, the disorder potential 
continues to be screened by the interactions. 

We thus have the paradoxical result that the disorder 
screening increases with increasing U, even though the 
overall current decreases. In contrast, the usual expecta- 
tion is that disorder screening should become worse as the 
Mott transition is approached, as occurs for strong disor- 
der. However, our result is consistent with earlier DMFT 
estimates that suggest that the elastic scattering rate due 
to weak disorder vanishes at the Mott transition. 26 Our 
result is also consistent with variational Monte Carlo cal- 
culations by Pezzoli et al.jfi 3 - who used a similar varia- 
tional approach to study the AHM. In particular, they 
introduced a variational renormalized site energy, similar 
to the parameter e used here, and found a monotonic re- 
duction of the effective disorder potential with increasing 
U. 

We note that the Gutzwiller calculations shown in 
Figs. S] and [S] break down above some critical interaction 
strength U C (W) at which the Mott MIT occurs. Where 
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FIG. 6: As above, but at quarter filling (n = 0.5). Results 
are for clean (W = 0) and disordered (W = D) systems for 
(a) a two dimensional lattice with L — 8 and D = 8t and (b) 
a one dimensional lattice with L — 50 and D — At. 



this happens depends on the level of disorder, and on 
the dimension of the lattice. For two dimensional clus- 
ters with W — 8t, variational calculations 23 find a first 
order Mott transition at U c ~ 12t, which is well into 
the regime where both strong correlations and disorder 
screening are significant. For one dimensional systems, 
the MIT happens at a smaller value; we have performed 
exact diagonalization calculations for small clusters with 
W — At that suggest that U c « 6t. It thus appears that 
the regime of strong screening and suppressed quasipar- 
ticle weight should be most easily observed in two and 
higher dimensions. 

For comparison, we show results for the persistent cur- 
rent at quarter filling in Fig. [SJ In this case, there is no 
Mott transition to worry about, and the Gutzwiller ap- 
proximation is expected to be meaningful over a large 
range of U values^ The most striking feature of this fig- 
ure is the diminished role of strong correlations at large 
U, meaning that J is not suppressed significantly at large 
U. 

As with the half-filled case, the figure shows that the 
PMGW current approaches the disorder-free current in 
the large U limit. The DFSGW current, on the other 
hand, is more weakly screened, and never approaches 
the clean limit value, but instead decays slowly from 
a maximum at U ~ W with increasing U. At quarter 
filling, then, the PMGW and DFSGW approximations 
are qualitatively different. In order to determine which 
result is correct, we have compared the two Gutzwiller 
results with published exact diagonalization results for 
the Drude weight at quarter filling in one dimension^ 
We find that the [/-dependence of the DFSGW current in 
Fig.[S[b) is surprisingly close to the exact diagonalization 
results, and in the remainder of this work we therefore 
show only DFSGW results at quarter filling. 

In summary, Figs. 0] and [5] illustrate one of the main 
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FIG. 7: Scaling of In J typ with system size in one dimension. 
Typical currents are shown for half filling with (a) U < W and 

(b) U > W. Data are for the PMGW approximation, solid 
lines are fits to the data from Eq. (1101) . Data are calculated 
for 10000/L disorder configurations, and error bars are give 
the statistical uncertainty in the typical current. Where error 
bars are not shown, they are smaller than the symbol size. 

(c) The scaling function for PMGW at half filling (left) and 
DFSGW at quarter filling (right). 



results of this paper, namely that within the Gutz wilier 
variational ansatz, strong correlations do not inhibit the 
screening of moderate disorder at half filling. 



B. Anderson Localization 

As discussed above, we can extract the localization 
length £ from the finite size scaling of J. In two di- 
mensions £ tends to be exponentially large 3 ^ for weak 
disorder, and we cannot reach the scaling limit L > £. 
We therefore illustrate our point by looking at the one 
dimensional case, where it is possible to obtain L»(, 
at least for U not too large. We showed in the previous 
section that the qualitative roles of screening and cor- 
relations are the same in one and two dimensions, and 
we therefore expect that the main lessons learned about 
localization in one dimension will also apply to higher 
dimensions. To explore the full range of behavior, from 
weak to strong correlations, we include results for large 
U where one dimensional systems are actually Mott in- 
sulating, but which are accessible in higher dimensions. 

Logarithmic plots of the typical current for PMGW 
calculations at half filling are shown in Fig. [Jj Lines are 
fits to the data of the form 



In J 



typ 



In J - - + — 



(10) 



where Jo, £, A, and y are fitting parameters. The last 
term is a finite size correction, meant to account for lead- 
ing order corrections when L is of order £. Clearly, the 



form of the correction cannot hold for small L, but it im- 
proves the quality of the fit significantly over the range 
of system sizes studied here. 

Two distinct trends can be seen in Fig. [7J for U < W, 
the magnitude of J typ increases with U, but the slope 
decreases; for U ^ W, the magnitude of Jt yp decreases 
with increasing U, and the slope continues to decrease. 
Since the slope is the inverse of the localization length, a 
decreasing slope corresponds to an increasing localization 
length. 

The localization length is shown in Fig. [8] At half 
filling, the PMGW and DFSGW results are qualitatively 
consistent; in both cases, the localization length grows 
monotonically with U. (Note that where £ exceeds the 
largest system sizes, namely L = 150, the scaling ansatz 
Eq. (|10p ceases to give quantitatively accurate values for 
£, and the results are qualitative only.) There is thus a 
region of the phase diagram where the overall magnitude 
of the current is suppressed by strong correlations, but 
the localization length grows with U. In contrast, at 
quarter filling, the localization length saturates at large 
U. 

In Fig. [7jc) , we plot the scaling function 



P 



din J 



typ 



rflnL 



(11) 



as a function of In J typ for both half filling and quarter 
filling. Note that the derivative in Eq. ([lip is obtained 
from the fitted curves in Figs. [7] (a) and (b). When one- 
parameter scaling holds, /? is a single- valued function of 
Jtyp. The main branch of the data shown in Fig. [TJc) 
corresponds to U < W and is, within the accuracy of 
our data, a function of J typ alone. There is also a small 
hook-shaped feature in the upper right hand corner of the 
figure that corresponds to a second branch which comes 
from U > W. This branch bends back towards the left, 
making (3 a multi- valued function of Jt yp - At half filling, 
it appears that one parameter scaling breaks down when 
U > W. In contrast, the quarter-filling data lies entirely 
on the main branch, and therefore appears to satisfy one- 
parameter scaling. 

This breakdown of one parameter scaling at half fill- 
ing follows directly from the two competing tendencies 
illustrated in Figs. 0] and [5] The main branch of the /3 
function comes from U < W, where J depends on U 
primarily through disorder screening, while the second 
branch comes from U > W where strong correlations 
are significant. While this breakdown may be hard to 
observe in one dimension because of the relatively low 
value of U c , it should be observable in two and higher 
dimensions. 

We note that our DFSGW results for £ at quarter fill- 
ing are entirely consistent with published DMRG calcu- 
lations by Nishimoto et ali££ It is difficult to compare 
results at half-filling, however, because the authors per- 
form their DMRG calculations for U > U C (W), and cal- 
culate £ as a function of filling, while enhanced disorder 
screening is predicted to occur at half-filling as U — >• U c . 
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FIG. 8: Localization length as a function of U for (a) PMGW 
and DFSGW at half-filling and (b) DFSGW at quarter-filling. 
PMHF results are shown for comparison. 



by an Aharonov-Bohm flux. 

We found that the persistent current is a nonmono- 
tonic function of U. For small U, the Coulomb interac- 
tions screen the disorder potential and the current is an 
increasing function of U. For large U, strong correlations 
suppress the current, which becomes a decreasing func- 
tion of U. The crossover between these two limits occurs 
at U » W. The main surprise in our calculations is that 
disorder screening persists when U is large. One conse- 
quence of this is that in the strongly correlated regime 
(U > W) the localization length is an increasing function 
of U, even though the current itself is a decreasing func- 
tion of U. We expect that this apparently-paradoxical be- 
havior should be observable two and higher dimensions. 



IV. CONCLUSIONS 



In summary, we have studied the effects of disorder and 
strong correlations on transport in the Anderson Hub- 
bard model. We have focused on the interaction-driven 
Mott transition at half filling with weak and intermedi- 
ate strength disorder. We used a Gutz wilier variational 
ansatz to describe the strongly correlated state near the 
Mott transition, and studied persistent currents induced 
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